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ABSTRACT 


In  a  multiple-input,  multiple-output  (M1M0)  radar  system,  two  or  more 
transmitters  emit  independent  waveforms,  with  the  resulting  refleetions  received  by  an 
array  of  receivers.  Recently,  MIMO  radar  has  become  a  subjeet  of  great  interest.  In  part, 
this  interest  is  due  to  the  potential  for  MIMO  techniques  to  reduce  radar  weight  and  cost, 
while  maintaining  performance  (as  compared  with  conventional  radar  approaches). 
However,  the  size  of  these  reductions  has  not  yet  been  quantified.  Likewise,  a  design 
process  that  minimizes  aperture  cost  (or  weight)  has  yet  to  be  developed. 

This  report  describes  a  process  for  designing  optimal  radar  apertures.  The  process 
treats  the  design  problem  as  one  of  minimizing  an  objective  function  under  performance 
constraints.  The  objective  function  is  based  upon  a  first-order  model  for  the  relationship 
between  cost  (or  weight)  and  performance,  and  is  derived  for  systems  employing  active, 
element-digitized  arrays.  A  systematic  process  for  optimizing  the  aperture’s  design  with 
respect  to  this  objective  function  is  presented,  and  equations  describing  the  optimal 
aperture  are  derived.  These  equations  provide  insight  into  the  optimal  relationship 
between  various  aperture  characteristics,  such  as  the  number  of  transmitters,  number  of 
receivers,  module  power  level,  and  virtual  array  length. 
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1.  INTRODUCTION 


In  recent  years,  Multiple-Input,  Multiple-Output  (MIMO)  radars  have  attracted  great  interest 
[1-15].  In  large  part,  this  interest  stems  from  the  potential  for  MIMO  techniques  to  yield  lighter-weight, 
lower-cost  airborne  and  surface-based  radar  systems  (as  compared  with  conventional  phased  arrays  of 
comparable  performance).  Other  benefits  include  increased  angle  [1-6]  and  Doppler  [1,7]  resolution, 
reduced  clutter  levels  [1,8],  sharper  moving  target  indication  (MTI)  clutter  notches  (yielding  lower 
minimum  detectable  velocity  [MDV])  [2,7],  lower  probability  of  intercept  [1,8],  and  relaxed  hardware 
requirements  [  1,8] — when  compared  with  conventional  radars  of  similar  complexity. 

These  benefits  all  stem  from  the  unique  fashion  in  which  MIMO  radars  transmit,  receive,  and 
process  signals.  In  MIMO  radar,  Nr  antennas  each  transmit  a  unique  waveform.  These  waveforms 
propagate  through  the  environment,  where  they  arc  reflected  by  targets  and  other  objects,  and 
subsequently  received  by  Nr  digital  receivers.  Faeh  receiver’s  output  is  then  processed  by  a  bank  of  Nr 
matched  filters  (i.e.,  one  matched  to  each  of  the  transmitted  waveforms).  As  such,  each  filter  essentially 
recovers  the  component  of  the  received  signal  due  to  a  single  transmitter.  Moreover,  the  phase  observed  at 
the  output  of  each  filter  behaves  like  that  of  a  signal  receiv  ed  by  virtual  sensor,  i.e.,  by  a  phantom  antenna 
element,  located  at  the  effective  phase  center  of  its  constituent  transmitter/reeeiver  pair.  Better  yet,  since 
there  are  a  total  of  Nt  ■  Nr  matched  filters  each  producing  data,  there  will  be  a  total  of  Nt  ■  Nf  output 
signals,  each  acting  like  it  originates  from  a  phantom  element  positioned  within  some  N \ ■  Nr  element 
“virtual  antenna  array.” 

As  with  any  conventional  phased  array,  the  N{  •  Nr  channels  that  form  the  MIMO  virtual  antenna 
array  can  be  combined  as  desired,  e.g.,  to  form  beams  or  maximize  signal-to-interference  plus  noise  ratio 
(SINR).  Not  surprisingly,  the  performance  achieved  by  combining  such  channels  will  depend  strongly 
upon  array  topology.  In  the  case  of  MIMO  radars,  in  particular,  it  is  the  virtual  array’s  topology  that  is 
paramount.  This  topology  depends  directly  upon  the  types  of  apertures  used  for  transmitting  and  receiving 
signals. 

To  date,  most  coherent  MIMO  radar  concepts  have  employed  one  of  two  apertures  types.  In  the 
first  type,  waveforms  are  transmitted  and  received  using  a  fully  filled  array,  as  shown  in  Figure  1(a)  (top). 
In  this  ease,  the  Nt  Ar  channels  form  a  virtual  array  that  is  densely  packed,  with  a  total  length  that  is 
nearly  twice  as  long  as  the  original  T/R  array.  This  type  of  v  irtual  array,  which  is  depicted  in  Figure  1(b) 
(top),  provides  the  flexibility  to  operate  in  both  MIMO  and  single-input  multiple-output  (SIMO)  modes, 
while  providing  enhanced  performance  in  MIMO  search  modes  along  with  excellent  MIMO  antenna 
pattern  control. 

In  contrast,  radar  applications  that  are  driven  by  aperture  length  and  spatial  resolution  (e.g., 
airborne  Ground  Moving  Target  Indicator  [GMTI]  radars)  often  prefer  a  second  type  of  MIMO  aperture. 
In  this  second  type  of  system,  waveforms  are  transmitted  and  received  using  a  pair  of  real  arrays,  see 
Figure  1(a)  (bottom).  One  of  these  arrays  is  large  and  highly  sparse  (i.e.,  undersampled  spatially),  while 


the  other  is  small  and  fully  filled  (i.e.,  Nyquist  sampled  spatially).  This  arrangement  is  attractive  heeause, 
with  appropriately  spaeed  elements,  the  resulting  virtual  array  will  be  both  very  long  and  fully  filled,  as 
shown  in  Figure  1(b)  (bottom).  Moreover,  sinee  the  virtual  aperture  is  uniformly  filled,  the  resulting 
antenna  patterns  will  be  well  behaved.  As  sueh,  MIMO  provides  sueh  systems  with  the  benefits  of  a  large 
aperture,  with  redueed  weight  and  eost,  and  without  the  elevated  sidelobe  levels  that  are  typieal  of  sparse 
arrays. 
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Figure  l.  (a)  Typical  MIMO  radar  apertures,  (b)  Associated  virtual  arrays.  Here,  "T/R  ”  denotes  an  integrated 
transmit/receive  element,  "Tx"  denotes  a  transmit-only  element,  and  "Rx  ”  denotes  a  receive-only  element. 


To  reduee  the  number  of  antennas  used  to  eover  a  desired  virtual  aperture,  this  seeond  type  of 
MIMO  aperture  is  preferred.  In  faet,  it  is  easy  to  show  that  the  number  of  antennas  will  be  minimized  by 
setting  the  number  of  transmit  and  reeeive  elements  to  be  equal.  As  sueh,  this  design  ehoiee  minimizes 
“complexity.”  Some  have  even  said  that  such  systems  will  “minimize”  eost  (or  weight).  However,  there  is 
no  simple  relationship  between  complexity  and  eost  (or  weight),  and  thus  no  quantitative  analysis  has 
been  performed  to  justify  this  conclusion.  Furthermore,  no  radar  system  design  procedure  has  been  given 
for  optimizing  the  design  of  sueh  arrays  with  respeet  to  eost  or  weight. 

To  address  this  problem,  we  outline  a  proeess  for  optimally  choosing  the  number  of  transmit  and 
reeeive  elements,  virtual  aperture  length,  and  transmitted  power  level,  under  constraints  on  seareh  rate 
and  detection  performance.  This  report  considers  both  conventional  (SIMO)  phased  array  radar  systems 
employing  filled  apertures  and  MIMO  radars  employing  sparse  apertures.  It  is  shown  that  the  “typieal” 
sparse  MIMO  design  approach  described  in  the  literature  (i.e.,  using  equal  number  of  transmit  and  reeeive 
elements)  is  usually  suboptimal.  Furthermore,  we  quantify  the  eosts  associated  with  optimal  MIMO  and 
SIMO  systems,  and  show  that  a  somewhat  modest  aperture  eost  savings  is  possible  by  using  an  optimally 
designed  MIMO  aperture,  with  larger  savings  achieved  when  sueh  systems  are  also  constrained  to  provide 
“exeess  spatial  resolution.”  These  eost  savings  will  be  offset,  to  some  degree,  by  the  increased  processing 
load  required  for  MIMO  systems.  Processing  eosts  are  not  addressed  in  this  report,  but  may  be  considered 
at  a  future  time. 
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2.  MIMO  RADAR  PRINCIPLES 


Consider  an  array  of  Nt  transmit  antennas  and  Nr  receive  antennas  (the  optimal  configuration  of 
these  antennas  is  discussed  in  Section  3),  In  MIMO  radars,  each  of  the  transmit  antennas  radiates  an 
independent  waveform.  Thus,  we  let  sn  denote  the  signal  radiated  by  transmit  antenna  n ,  where 
A7  =  l,...,jV,  (he.,  sw  is  a  Pwx  1  complex  vector  containing  baseband  samples  of  the  /?th  radiated  signal). 
Furthermore,  let  0t  represent  various  parameters  (e.g.,  location  coordinates)  associated  with  target  /,  and 
let  a rv(^)  be  the  corresponding  Nf  x  1  response  vector  for  the  transmit  antenna  array.  Under  these 
conditions,  the  signal  reflected  from  point-target  /  is  proportional  to  aA  (6>  )S,  where  the  jV,  x/1  signal 
matrix  S  is  defined  via  S=[s,  s2  •  -  s,v  J  . 

The  reflected  signal  propagates  through  the  environment  and  is  received  by  an  array  of  antennas. 
Assuming  each  receiver  collects  R  »  Pw  samples,  the  Nr  x  R  sample  matrix  at  the  output  of  the 
receive  array  has  the  form 

X  =  Za(^)a«<(^)ar.(^)S;  +  Z“(^)a«-(^)ar,(^)Sc  +:N  -  U) 

/  v 

where  subscripts  t  and  c  index  into  the  set  of  targets  and  clutter  patches,  respectively.  Here,  a ( 9 ( ) 
represents  a  complex  scale  factor  associated  with  target  /,  while  a /tr(^/)  denotes  the  receive  array's 
yVr  x  1  response  vector  for  target  t.  The  Nt  x  R  matrix  S;,  on  the  other  hand,  describes  how  the  various 
waveforms  propagate  to,  and  are  reflected  by,  target  /.  That  is,  the  nh  column  of  S,  contains  the 
convolution  of  sn  with  the  target's  impulse  response.  Furthermore,  since  R  »  PH  ,  Sf  will  also  contain 
many  zeros  to  account  for  propagation  delays  to/from  the  target.  Variables  a{6  ),  (#  ) ,  aA  (tf  ) ,  and 

S  are  defined  similarly  with  respect  to  each  clutter  patch.  Lastly,  the  Nr  x  R  matrix  N  represents  noise 
as  well  as  other  interference. 

After  being  received,  the  sample  matrix  is  processed  to  extract  target  detections  and/or  to  form 
images.  Target  detection,  for  example,  can  be  performed  by  passing  the  received  samples  through  a  bank 
of  matched  filters.  Often,  this  filtering  is  implemented  in  stages,  wherein  the  initial  stage  applies  a  set  of 
"‘waveform’'  matched  filters  to  form  N{  x  R'  matrices, 

V  = 

m 


where  \'m  denotes  the  m  1  row  of  X,  *  denotes  convolution,  and  h(i  denotes  the  matched  filter  for  signal 
sn .  After  matched  filtering,  the  rows  of  the  N,Nr  x  R'  matrix  V’  =  [~ V'/  •••  Ysr  J  represent  N,Nr 

channels  of  data.  It  is  thus  possible  to  linearly  combine  these  channels. 


*h r 


*h^, 


,  m  =  l,...,N 


(2) 


3 


z  =  w,,Y  » 


(3) 


to  achieve  some  desired  goal,  such  as  the  formation  of  beams  [1,3,8]  or  the  maximization  of  SINR  [2,3]. 
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3.  SEARCH  RADAR  APERTURE  OPTIMIZATION 


In  this  section,  the  problem  of  designing  optimal  search  radar  apertures  is  addressed.  Optimization 
is  performed  from  a  cost  perspective,  under  performance  constraints.  Comparisons  between  optimized 
MIMO  and  SIMO  systems  are  presented. 

3.1  OPTIMIZED  MIMO  SEARCH  RADAR 

To  begin,  let  us  assume  the  use  of  a  critically  sampled  MIMO  virtual  aperture  of  length  L 
Assuming  each  element  has  physical  aperture  length  of  A/2 ,  Nt  is  given  by 


Nt  =2L/(A-A'r)  .  (4) 

Next,  impose  constraints  upon  the  search  performance  of  this  system  by  requiring  it  to  have  a  predefined 
power-aperture  (PA)  product.  For  MIMO  systems,  the  PA  product  is  proportional  to£  ,  where 

f  =  pm.  Nt.Nr  ,  (5) 

and  Pm  is  the  average  RF  power  radiated  from  each  transmit  antenna.  Substituting  our  expression  for  Nt 
into  the  PA  product,  we  have 

£  =  2 -L-PJA.  (6) 

Assuming  the  length  of  the  MIMO  virtual  aperture  is  a  free  variable  (i.e.,  the  resolution  is  unconstrained), 
we  can  solve  for  L  as  a  function  of  £  ,  yielding 

L  =  &/{2Pm).  (7) 

Our  goal,  then,  is  to  choose  AT. ,  Nt,  and  Pm  so  as  to  minimize  radar  aperture  cost  under  the 
performance  constraint  above.  To  accomplish  this  goal,  we  first  define  a  model  relating  aperture  design 
characteristics  to  cost.  Table  1  lists  various  aperture  components,  and  describes  their  contribution  to  cost. 
In  particular,  each  row  in  Table  1  describes  a  specific  component  of  the  system.  Within  the  row,  columns 
describe  how  the  component’s  cost  is  related  to  fundamental  aperture  characteristics.  For  example, 
consider  the  antenna  receiver  (Rx)  module,  which  includes  the  low  noise  amplifier  (LNA)  and  other 
devices.  As  Nr  increases,  the  number  of  such  modules  must  also  increase  (for  both  MIMO  and  SIMO 
systems),  thereby  raising  their  contribution  to  overall  aperture  cost  proportionally.  In  Table  1,  this 
relationship  is  shown  by  the  entering  “  Nr  ”  in  the  columns  labeled  “Dependence  on  Nt  or  Nr  '' 

Similarly,  other  component  costs  (c.g.,  the  low  power  transmit  circuitry)  will  depend  solely  on  a 
single  aperture  characteristic  such  as  Nt  or  NtPm  (i.e.,  total  radiated  power).  However,  not  all  component 
costs  will  depend  on  a  single  characteristic.  Some  costs  will  depend  on  multiple  factors,  for  example,  the 
cost  associated  with  the  high  power  portion  of  each  transmit  module  (which  includes  the  high-power 
amplifier  and  other  circuitry)  will  add  to  the  total  aperture  cost  in  a  way  that  depends  both  on  N(  and 
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Pm/Em  (i.e.,  the  power  to  each  module).  Likewise,  the  cost  to  cool  each  module  will  depend  on  three 
factors:  Nt,  (l  -Em)Pm/Em  (i.e.,  the  power  dissipated  by  the  module),  and  (l  -  Em)Pm/(Em  -  L  -  H)  (i.e., 
the  power  density  of  the  array,  in  watts  dissipated  per  square  meter,  assuming  an  Lx  II  tile-based  array). 

For  the  sake  of  simplifying  the  analysis  that  follows,  we  assume  (1)  the  use  of  a  DC-to-DC 
converter  at  each  transmit  element,  and  (2)  the  exponent  values  ( a ,  /?,  c,  and  <J)  in  Table  1  are  equal  to  1. 
These  assumptions  lead  us  to  a  closed  form  solution  for  the  optimum  array  design,  avoiding  the  need  for 
numerical  solutions.  In  practice,  though,  a  more  accurate  model  can  be  obtained  by  removing  these 
assumptions. 


Table  1 

Cost  Model  for  SIMO  and  MIMO  Radar  Apertures* 


SUBSYSTEM 

COMPONENT 

% 

Fixed 

Term 

Dependence  on  N  t  or  N r 

SIMO  Model  MIMO  Model 

Dependence  on 
Power 
(times  Nf) 

Dependence  on 
Power  Density 
(times  N,) 

Antenna 

Radiator 

Nt 

Nt,Nr 

- 

Structure 

Yes 

N. 

(Nt,Nr) 

- 

High  Power Tx 

Nt 

N, 

(PmIEm)9 

Low  Power  Tx 

Nt 

Nt 

- 

Rx 

N, 

Nr 

- 

DC-DC  Conversion 

Variable 

Variable 

Module  Cooling 

Nt 

Nt 

(i  -EmypjEm 

(1  -Em)Pml 
{Em*L'H) 

Receiver/ 

Exciter 

Centralized  Exciter 

Yes 

- 

- 

- 

- 

Distributed  Exciter 

N  t 

Nt 

- 

- 

Distributed  Receiver 

Nr 

Nr 

- 

Power  & 
Cooling 

Power  Conditioning 

- 

- 

(PmlEsf 

- 

Chilling 

- 

- 

(Pm(1-£s)/£s)° 

- 

Overhead 

Test  &  Evaluation 

Yes 

- 

- 

- 

- 

Management 

Yes 

" 

" 

" 

- 

*Where  Em  denotes  module  efficiency  and  Es  is  system  efficiency.  The  model  assumes:  (1)  active,  tile-based 
transmit  arrays  and  element-digitized  receive  arrays,  (2)  available  DC  power,  (3)  noncontiguous  sparse 
antenna  structure,  and  (4)  module  cooling  costs  that  depend  on  power  and/or  power  density. 


Continuing  under  these  assumptions,  Table  1  suggests  the  following  parametric  model  for  aperture 
hardware  cost: 


(c,  +  CpPm  +  CdPjL)N,  +  CrNr  +  Cf  .  (8) 

The  total  cost  to  build  and  test  the  aperture,  in  contrast,  is  modeled  as  the  sum  of  the  array  hardware  cost 
(8)  plus  a  percentage  to  cover  test,  evaluation,  and  management.  Nonetheless,  the  best  set  of  design 
parameters  will  be  those  that  minimize  (8). 

Note  that  the  coefficients  in  our  aperture  cost  model  (8)  play  an  important  role:  they  quantify 
relationships  between  various  aperture  design  characteristics  and  aperture  cost.  The  coefficient  Cf ,  for 
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example,  represents  the  total  fixed  cost  (i.e.,  the  portion  of  eost  that  is  independent  of  the  aperture  size, 
power  level,  or  power  density).  Likewise,  the  coefficients  Cr,  C, ,  Cp ,  and  Ct/  arc  used  to  describe  the 
dependence  of  eost  on  Nr,  Nt ,  and  N(Pm/L ,  respectively  (where  L  is  the  length  of  the  active 

portion  of  the  aperture;  array  height  H  is  assumed  to  be  fixed).  The  actual  values  assigned  to  these 
coefficients  (i.e.,  Cf  ,  Cr ,  Ct ,  C  ,  and  C, )  will  be,  in  practice,  determined  empirically  (e.g.,  by  itemizing 
the  cost  for  each  component  in  an  existing  radar  system,  then  factoring  these  costs  into  the  form  given  by 
(8)). 

Furthermore,  optimizing  with  respect  to  the  parametric  form  in  (8)  allows  us  to  investigate  certain 
special  eases  directly.  For  example,  we  can  set  Ca  =  0  in  (8),  yielding  a  simpler  eost  model.  Setting  Cti 
0  implies  the  cost  to  cool  the  aperture  does  not  depend  on  power  density.  In  contrast,  allowing  Ci{  to  be 
nonzero  allows  one  to  investigate,  to  a  degree,  the  impact  of  higher  eost  cooling  technologies  (e.g.,  liquid 
cooling  versus  forced  air)  needed  for  high  power  density  radars.  Note  that  (8)  models  the  relationship 
between  eost  and  power  density  as  linear.  In  practice,  this  is  an  approximation  as  illustrated  in  Figure  2. 
Additional  accuracy  is  possible  by  allowing  costs  associated  with  cooling  to  grow  nonlinearly  with  power 
density.  Such  models,  however,  typically  require  numerical  solution. 


</> 

o 

o 


Power  Density 


Figure  2.  First  order  model  relating  cost  to  power  density/cool  if  ig  technology 


Substituting  our  expressions  for  i\\  and  L  above,  we  may  rewrite  our  eost  funetion  dircetly  in  terms 


of  N  and  P_  : 


r 


Cost  {Nr,Pm)  =  Cf+  Cr  Nr  + 


(l» 


The  partial  derivatives  with  respeet  to  Nr  and  Pm  are 


(10) 
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and 


J-CnitlN  P\  4C^+C^A  2Ct,Pj  +  Cp&Pm  +  C,& 
dP  \  r'  m>  N  PA  N.PlA 


(11) 


To  find  the  minimum  cost  system,  we  set  these  partial  derivatives  equal  to  zero  and  solve  for  Nr  and  Pm 
as  a  function  of  Ct ,  Cp ,  Cd ,  and  A .  One  solution  is 


and 


(12) 


(13) 


Note  that  there  are  several  other  solutions.  Consequently,  we  need  to  verify  that  the  values  above  produce 
a  cost  minima  by  computing  the  second  derivatives: 

A  52  ,  ■.  4C.P2  +  2C  £AP  +  2C,£A 

f  t  ——r Cost (Nr,P  )  =  -  — - £__• - — 

SN 2  V  ’  N^PA 


fyy=  C0St{N,,Pm)  = 


4CyS+2C£U>„+2C,#  »C„P.  +  2CM  4  C, 


Nrpy. 


NrP, U 


NP  A 


f  A 

J  xy  — 


_5 _ d_ 

dp'dN 


Cost(Nr,Pm)  = 


2  CdPi  +  Cp&Pm  +  C,a  4CjPm  +  C,£l 


n]p:A 


n2p  a 


(14) 


Clearly,  >  0  .  From  calculus,  then,  our  solutions  for  Nr  and  Pm  will  provide  minimum  cost  as  long  as 
fxxfy  -  fl  is  greater  than  zero.  The  second  derivative  test  f>y  -  f*,  simplifies  to 


8CjCf2^(32C,C,2+(C^  +  24QC/,C,  )pC,ZA/C,  +  l2C2pC,q) 
C,(^C^A+C^ 


(15) 
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Moreover,  since  C, ,  Cr,  Cp ,  Q,  /i,  and  £  are  all  positive  real  numbers,  fxxfyy  -  fxv  must  be  greater 
than  zero.  Hence,  our  equations  for  Nr  and  Pm  describe  an  optimum  cost  MEMO  search  radar. 

Next,  the  cost  ot  this  system  can  be  found  by  substituting  our  solutions  back  into  the  cost  function. 
This  cost  is 


(16) 


Finally,  the  optimal  values  for  N{  and  L  can  be  found  by  substituting  our  solutions  for  Nr  and  Pm  into 
our  original  constraint  equations.  The  optimal  values  for  Nt  and  L  are 


N,  - 


S 


and 


U8CCJA +c,a 

lea 

V  c ^  V 

2  Cd 

L  = 


& 


pC'ZA/C, 


(17) 


(18) 


Note  that  the  optimal  ratio  of  Nr  to  Nr  is 


•K-.  +  c.Va^i-c, 

2Cr  ' 

What  can  be  said  about  this  ratio?  Well,  let's  start  by  considering  what  happens  when  Ct  ~  Cr . 
(This  ease  is  a  reasonable  starting  point  because  many  of  the  costs  that  depend  on  Nt  and  Nr ,  such  as 
radiators  and  structures,  will  be  identical,  while  other  components  used  to  build  Rx-only  channels  will  be 
of  similar  cost  to  those  used  in  building  the  low  power  parts  of  the  Tx-channels).  In  this  special  case,  our 
analysis  gives  us  an  important  insight: 

With  the  chosen  cost  model  and  constraints ,  the  number  of  receive  elements  in  an  optimal  MI  MO  search 
radar  should  he  at  least  double  the  number  of  transmit  elements. 
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These  equations  also  provide  another  bit  of  insight.  Consider  the  limit  as  costs  associated  with 
“power  density”  go  to  zero  (i.e.,  Cd  goes  to  0).  In  this  limit,  it  doesn’t  cost  any  more  to  cool  a  small  array 
than  it  does  to  cool  a  large  one.  Moreover,  the  optimal  number  of  receivers  will  approach 


N_  = 


(20) 


Likewise,  the  number  of  transmitters  will  go  to  one  (it  cannot  go  below  one!).  Thus,  as  Cd  goes  to  zero, 
the  best  MIMO  radar  is  a  “ubiquitous  radar” — to  borrow  a  term  popularized  by  Merrill  Skolnik  in  [16]. 
Hence,  MIMO  radar  can  be  viewed  as  a  generalization  of  the  ubiquitous  radar  concept,  allowing  active 
array  antennas  with  finite  cooling  costs  to  provide  the  benefits  of  “ubiquitous  radar”  operation  [8]. 

Before  moving  to  the  next  section,  some  final  notes.  First,  the  equations  above  do  not  explicitly 
model  antenna  elevation  pattern  gain.  However,  one  should  be  able  to  simply  adjust  the  PA  product  to 
allow  for  the  elevation  gain.  Second,  the  above  analysis  is  approximate  in  the  sense  that  it  assumes  the 
dissipated  power  density  is  uniform  over  the  full  virtual  aperture.  For  MEMO  systems  employing  a  large, 
sparse  transmit  array  (and  a  small  filled  receive  array),  this  is  approximately  true.  However,  for  systems 
employing  a  small  filled  transmit  array  and  a  sparse  receive  array,  it  is  unlikely  unless  the  module  power 
is  very'  low.  However,  under  the  cost  model  above,  the  use  of  small  filled  transmit  arrays  will  be  more 
costly  than  sparse  transmit  arrays.  Since  we  seek  the  minimum  cost  MIMO  system,  we  need  only 
consider  the  sparse  transmit  array  case. 

3.2  OPTIMIZED  SEARCH  RADAR  BASED  ON  CONVENTIONAL  BEAMFORMING 

Next,  let’s  compare  the  optimized  MIMO  search  radar  (above)  to  a  conventional  radar  system  that 
has  been  optimized  for  the  same  purpose.  In  the  ease  of  the  conventional  SIMO  radar,  a  single  aperture  of 

length  L  is  assumed  to  be  shared  for  transmit  and  receive.  Consequently,  Nt  =  Nr  w  ith 

Nt=  2L/A  +  1  Nr  =  2L/X  + 1  .  (21) 

As  above,  we  constrain  search  performance  via  Pm  •  Nt  •  Nr  =  £  (or  equivalently  Pm  -N2  =  £ ).  Then, 

substituting  our  expression  for  Nt  into  the  PA  product,  we  have 

L  =  {1/ 2)(V?7^-l)  •  (22) 


Assuming  the  cost  function  has  the  same  form,  we  substitute  our  expressions  for  Nn  Nr,  and  L  to  get 


r,p >  c,fs+c,f^c,4F„+crP„fs  icdp„4s 


(23) 
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and 


cpj4  cf  cr^  +  c,^  +  cfyl^  +  cppm^ 


2 cJl  ,  cdJFmJj 


(24) 


To  solve  for  the  minima  in  closed  form,  we  need  to  equate  this  derivative  to  zero  and  solve  for  Pm , 
as  we  did  for  the  MIMO  case.  In  practice,  this  process  is  tricky.  To  approach  this  problem,  we  could 
collect  the  Pm  terms  of  dCost[Pm)jdPm  .  This  gives  a  fraction  and,  since  we  are  setting  this  fraction  equal 
to  zero,  wc  can  discard  the  denominator.  The  remaining  term  (i.e.,  numerator)  is  a  generalized  polynomial 
having  components  involving  Pm  raised  to  various  powers  (i.e.,  2,  1.5,  1,  and  1/2).  So,  we  would  next  do 
a  change  of  variables,  letting  x  =  y[P^ .  This  gives  us  a  quartic  equation  (i.e.,  x  is  raised  to  the  1st,  2nd,  3rd, 
and  4th  powers).  The  closed-form  solution  to  such  quailic  equations,  although  known,  is  too  long  to 
provide  here  (Mathcmatica™  can  provide  the  solution).  As  a  reasonable  alternative,  one  can  simply  solve 
numerically  for  the  roots.  This  is  the  approach  used  in  Section  3.3  below. 

3.3  SEARCH  RADAR  DESIGN  EXAMPLE 

To  illustrate  the  optimization  approaches  described  in  Sections  3.1  and  3.2,  wc  require  two  sets  of 
input  parameters,  (1)  a  radar  system  design  specification  (as  defined  by  £  and  A),  and  (2)  a  set  of  cost 
model  coefficients  (C,  ,Cr ,  Cp,  C,,  and  Cf  ).  As  a  hypothetical  example,  wc  assume  -  le4,  A  =  .02, 
C,  =  Cr  =  7200,  Cp  =  732,  Cd  =  1200/. 2  (aperture  height  .2  m),  and  C f  =  2.8e6.  In  addition,  8%  of  the 
budget  is  reserved  for  test,  evaluation,  and  management.  Note  that  these  parameters  are  meant  only  to 
loosely  represent  costs  for  a  low-production-volume  microwave  radar  system. 

Inserting  these  parameters  into  our  MIMO  model  (Section  3.1),  the  optimal  number  of  receive 
elements  (12)  is  Nr  =  53.314.  Likewise,  the  module  power  (13)  is  Pm  =  10.954,  the  virtual  aperture 
length  (18)  is  /  =  9.129,  and  the  number  of  transmit  elements  (17)  is  Nt  =  17.122.  Consequently,  the 
expected  cost  of  the  array  hardware,  as  defined  by  (16),  will  be  S3.57M  excluding  overhead.  Likewise, 
the  cost  including  overhead  will  be  S3.88M. 

Of  course,  we  may  insert  these  same  parameters  and  coefficients  into  the  S1MO  model  of  Section 
3.2.  In  doing  so,  we  find  that  the  derivative  (24)  equals  zero  at  Pm  —  1.076.  From  (21)  (23),  it  follows 
that  the  optimal  S1MO  radar  is  characterized  by  L  =  0.954,  Nr  =  96.387,  and  an  array  cost  (including 
overhead)  of  S5.34M. 

Two  observations: 

1 .  The  MIMO  aperture  costs  less  than  the  SIMO  aperture. 

2.  The  comparison  is  potentially  flawed  because  it  assumes  the  same  cost  coefficient  values  for 


MIMO  and  SIMO.  This  assumption  is  reasonable  only  if  the  SIMO  radar  uses  separate  transmit 
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and  receive  arrays.  In  this  case,  however,  there  is  no  reason  why  the  transmit  and  receive 
apertures  would  be  equal  size.  In  faet,  if  separate  transmit  and  reeeive  arrays  arc  used  in  a 
SIMO  system,  many  of  the  same  cost  advantages  will  be  realized  (note:  the  SIMO  system 
would  still  have  a  higher  power  density  and  would  require  more  eooling — whieh  is  potentially  a 
significant  eost  faetor). 

Observation  2  ean  be  addressed  via  modification  of  the  SIMO  eost  model.  Today,  most 
conventional  SIMO  radars  use  the  same  antennas  used  for  Tx  and  Rx.  In  this  case,  a  single  module  can  be 
used  at  each  antenna  location.  This  module  will  be  more  sophisticated  than  the  MIMO  transmit  (or 
reeeive)  modules  assumed  above.  For  example,  it  will  require  a  T/R  switeh  or  eireulator,  along  with 
transmit  phase  and  amplitude  eontrol.  These  additional  components  will  increase  individual  module  cost, 
but  this  increase  will  be  offset  by  the  need  for  a  single  T/R  module  rather  than  two  separate  modules  (i.e., 
Tx-only  and  Rx-only).  In  faet,  assuming  a  single  T/R  module  is  used  instead  of  two  simpler  modules,  the 
SIMO  eost  function  will  be 


Cost{Pm)  =  Cf 


,  c,Js+cppmJl 


1£AJL 


(25) 


Furthermore,  allowing  the  T/R  module  eost  parameter  to  rise  to  Ct  =  8900,  the  partial  derivative  of 
Equation  (25)  will  equal  zero  at  Pm  =  0.777  W.  From  Equations  (21),  (22),  and  (25),  this  implies  the 
optimal  SIMO  radar  will  be  eharaeterized  by  L  =  1.124,  Nr  =  113.446,  and  a  resulting  array  eost, 
including  overhead,  of  S4.72M. 

Note  that,  even  with  the  sharing  of  transmit  and  receive  hardware,  conventional  radar  eosts  arc 
expected  to  significantly  exceed  MIMO  radar  costs. 
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4.  APERTURE  OPTIMIZATION  WITH  RESOLUTION  CONSTRAINT 


In  this  section,  we  address  the  problem  of  designing  optimal  search  radars  providing  a  specified 
spatial  resolution.  As  in  the  previous  section,  we  consider  both  MIMO  and  SIMO  systems. 

4.1  OPTIMIZED  MIMO  SEARCH  RADAR  WITH  RESOLUTION  CONSTRAINT 

First,  we  consider  the  design  of  MIMO  radars  that  minimize  cost  under  both  power-aperture  and 
resolution  constraints.  As  in  our  earlier  MIMO  discussion,  a  critically  sampled  MIMO  virtual  array  is 
assumed.  Hence,  Nt,  Pm ,  and  the  cost  function  are  as  defined  in  Section  3.1.  This  time,  however,  we 
begin  by  substituting  our  expressions  for  Nt  and  Pm  into  our  cost  function,  yielding 


Cost(Nr)  =  CrNr  +  Cj 


+^A+^A+1(^L 

Nr  NrL  NrA 


(26) 


Next,  since  resolution  is  constrained  (i.e.,  L  is  given),  we  simply  choose  Nr  to  minimize 
cost.  The  derivative  of  our  cost  function  with  respect  to  Nr  is 


8 

<3/Vr 


Cost(Nr)  =  Cr 


2CJ:  +  C p8 A  L  +  Cd8A 
N;AL 


(27) 


The  derivative  equals  zero  when 


12  C,L2  +  Ct£AL  +  C^A 

CrLA 

Likewise,  the  second  derivative  of  our  cost  function  is 


(28) 


fM 


Cost  ( Nr ) 


2CJ  +  2CpLZ  |  4 C,L 
NlL  NlA 


(29) 


At  the  critical  point  of  Equation  (28),  it  can  be  shown  that  f„(Nr )  is  positive.  Thus,  the  critical  point 
gives  the  value  of  Nr  associated  with  the  minimum  cost  design  solution.  Likewise,  substituting  Equation 
(28)  into  Equation  (4)  yields  the  optimal  Nt,  with  Pm  computed  directly  from  Equation  (6). 


Note  that  the  optimal  ratio  Nj N,  is  now  given  by 


C,  f  C^A  +  CpUA 
Cr+~  2  CrL2 
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Furthermore,  this  ratio  is  directly  related  to  the  optimal  unconstrained  MIMO  system  through  Equation 
(18),  i.e.,  substituting  Equation  (18)  into  Equation  (30)  yields  Equation  (19). 

Lastly,  observe  that  as  L  grows  larger  in  Equation  (30),  the  eosts  associated  with  eooling  beeomc 
less  significant.  Likewise,  the  module  power  decreases  (to  a  point  which  it  is  dominated  by  the  load  from 
the  low-power  electronics),  reducing  costs  associated  with  high-power  amplification,  DC  conversion, 
power  conditioning,  and  chilling.  In  the  limit,  as  L  grows  very  large,  the  Ct  and  Cr  cost  terms  dominate 
the  optimization.  Furthermore,  if  C,  ~  Cr ,  then  (in  the  limit)  the  optimum  system  will  have  Nr  ^  Nr .  In 
this  sense,  the  system  will  be  “balanced.”  However,  the  aperture  size  required  to  approach  this  limit  can 
be  quite  large,  making  this  limiting  scenario  the  exception,  not  the  rule. 

4.2  OPTIMIZED  SIMO  SEARCH  RADAR  WITH  RESOLUTION  CONSTRAINT 

Finally,  let’s  eompare  our  optimized  MEMO  search  radar,  with  resolution  constraints,  to  a 
conventional  SIMO  system  designed  for  the  same  purpose.  Here,  Nt  and  Nr  are  defined  by  Equation 

(21)  and  the  power-aperture  product  is  constrained  via  Equation  (5).  With  the  resolution  also  constrained 
(i.e.,  L  given),  the  only  unknown  is  Pm  .  Thus,  by  substituting  Equation  (21)  into  Equation  (5),  we  have 

Pm=&2/(2L  +  A)2  .  (31) 

The  cost  is  then  computed  directly  from  Equation  (25)  or  (23),  depending  upon  whether  a  single  aperture 
(i.e.,  shared  transmit  and  receive  array)  or  pair  of  apertures  is  used. 

4.3  RESOLUTION  CONSTRAINED  SEARCH  RADAR  EXAMPLE 

To  illustrate  the  incorporation  of  resolution  constraints,  let’s  assume  the  same  values  for  C,,  Cp,  Cj, 
and  Cr  as  in  Section  3.1.  Furthermore,  suppose  we  constrain  L  =  25  meters  (the  approximate  useable 
length  of  a  narrow-body  airplane,  such  as  737,  A-318,  or  El  70).  In  this  case,  Equation  (28)  dictates  that 
the  optimum  MIMO  radar  system  will  have  Nr  =  62.048.  This  implies  N{  =  40.291,  Pm  =  4,  and  the 
projected  array  cost  (including  overhead)  for  this  system  will  be  S4.02M. 

Likewise,  assuming  the  same  values  for  C, ,  C  ,  Cd ,  Cr,  and  L,  Equation  (21)  dictates  the 
optimum  SIMO  radar  will  have  Nt  =  Nr  =  2501.  Furthermore,  Equation  (31)  dictates  Pm  =  1.599e-3, 
which  results  in  a  projected  array  eost  (ineluding  overhead)  of  S42.19M. 

Of  course,  as  discussed  in  Section  3.2,  it  is  common  to  implement  SIMO  radars  using  a  shared  T/R 
aperture  and  combined  T/R  modules  (rather  than  separate  Tx  and  Rx  modules).  In  this  case,  the  optimum 
SIMO  radar  will  still  have  Nt  =  Nr  =  2501  and  Pm  =  1.599c-3.  Moreover,  the  projected  array  cost 
(including  overhead)  drops  to  S27.24M  because  of  efficiencies  gained  by  using  combined  T/R  modules 
and  a  shared  T/R  aperture.  However,  the  SIMO  radar  cost  is  still  nearly  an  order  of  magnitude  greater 
than  the  MIMO  system  cost,  while  providing  the  same  resolution  and  PA  produet.  This  higher  eost  is 
partially  because  we  required  an  element  digital  SIMO  array  architecture.  Today,  large  SIMO  radars 
usually  employ  analog  subarray  beamforming.  Using  analog  subarrays  will  further  reduee  the  SIMO 
radar’s  cost. 
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5.  WEIGHT  OPTIMIZATION 


In  Sections  3  and  4,  our  analysis  focused  explicitly  on  minimizing  aperture  cost.  Clearly,  cost 
minimization  is  not  always  the  designer's  main  objective.  For  example,  systems  that  will  fly  on  small 
aircraft  may  be  driven  by  the  need  to  minimize  payload  weight  rather  than  cost.  Nonetheless,  the  first 
order  objective  function  will  still  take  the  form  given  by  Equation  (8).  Hence,  even  though  the  objective 
function  coefficient  values  will  change,  our  optimization  approach,  as  well  as  the  equations  describing  the 
solution,  should  still  apply. 
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6.  SUMMARY  AND  CONCLUSIONS 


Applying  traditional  engineering  techniques,  it  has  been  shown  that  it  is  possible  to  optimize 
MIMO  and  S1MO  apertures  to  minimize  an  objective  (e,g,,  eost  or  weight  function)  subjeet  to  constraints 
on  performance.  A  simple  procedure  was  presented.  Although  simplistic,  our  procedure  provided  useful 
insight  into  the  amount  of  savings  that  ean  be  achieved  by  using  MIMO  technology.  It  also  provided 
insights  into  various  MIMO  radar  design  relationships.  For  example,  under  the  assumed  hardware  and 
cost  model  (for  aetive  element-digital  arrays),  MIMO  search  apertures  typieally  have  Nr  »  Nt. 

To  illustrate  our  approach,  wc  presented  design  examples  assuming  £  =  10000,  as  well  as  various 
hypothetical  eost  parameters  (meant  to  represent  a  low-produetion  mierowave  radar).  Table  2  summarizes 
the  resulting  system  designs. 


Table  2 

Summary  of  SIMO  and  MIMO  Radar  Design  Examples 


NO.  OF 

NO.  OF 

LENGTH 

MODULE 

APERTURE 

CONSTRAINTS 

TRANSMITTERS 

RECEIVERS 

(REAL/VIRTUAL) 

POWER 

COST  ($) 

cinn 

PA 

96 

96 

0.95  m 

1.1  Watts 

$5.3M 

ol  IV1U 

113 

113 

1.1  m 

.8  Watts 

S4.7M 

MIMO 

PA 

17 

53 

9.1  m 

10,9  Watts 

$3  9M 

SIMO 

PA  &  Resolution 

2501 

2501 

25  m 

1 .6  mWatts 

S42.2M 

2501 

2501 

25  m 

1 .6  mWatts 

527.2M 

MIMO 

PA  &  Resolution 

42 

62 

25  m 

4  Watts 

$4.0M 

(1)  Separate  T  and  R  modules  and  antennas 

(2)  Combined  T/R  module  and  antenna 


Our  example  suggests  that  a  power-aperture  (PA)  constrained  MIMO  aperture  will  be  moderately 
less  expensive  as  compared  with  a  comparable  SIMO  aperture.  In  contrast,  when  both  power-aperture  and 
resolution  are  constrained,  MIMO  apertures  ean  be  dramatically  less  costly.  It  should  be  noted,  however, 
that  MIMO  systems  require  additional  digital  processing  (not  considered  here),  MIMO  radars  may  also 
saerifiee  some  flexibility,  which  is  important  in  applications  requiring  multimode  or  multifunction 
operation.  Furthermore,  “critically  sampled”  MIMO  radars  will  be  more  degraded  by  module  failures  than 
comparable  SIMO  active  electronically  scanned  array  (AESA)  radars.  As  a  result,  additional  eost  may  be 
incurred  to  incorporate  fault  tolerance.  Finally,  limited  waveform  orthogonality  is  known  to  degrade 
MIMO  radar  system  performance  in  a  way  that  is  unique  to  MIMO  radars.  Hence,  eost  is  only  one  factor 
that  should  be  considered  in  selection  of  a  radar  architecture. 
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